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We propose a model for planar nanoscale dielectric capacitor consisting of a single layer, insulating 
hexagonal boron nitride (BN) stripe placed between two metallic graphene stripes, all forming 
commensurately a single atomic plane. First-principles density functional calculations on these 
nanoscale capacitors for different levels of charging and different widths of graphene - BN stripes 
mark high gravimetric capacitance values, which are comparable to those of supercapacitors made 
from other carbon based materials. Present nanocapacitor model allows the fabrication of series, 
parallel and mixed combinations which offer potential applications in 2D flexible nanoelectronics, 
energy storage and heat-pressure sensing systems. 

PACS numbers: 


I. INTRODUCTION 

While sustainable clean energy production and stor¬ 
age have been one of the critical concerns of our century, 
nanoscale capacitors are developing as an energy storage 
medium due to their recyclable high-energy storage based 
on charge-separation.® In particular, their flexibility, light 
weight and high performance have made nanoscale ca¬ 
pacitors advantageous in the field of electronics and en¬ 
ergy storage. With this regard, graphene with its one- 
atom thick layer and high chemical stabilitjl^ has been 
recently used in electrochemical capacitors.®^!® As op¬ 
posed to electrochemical double-layer capacitors or su¬ 
percapacitors, where an aqueous solution is used to store 
electric energy with electron transfer through chemical 
reactions, a nanoscale dielectric capacitor (NDC) uses a 
solid dielectric between the metallic plates and can attain 
properties that are superior to other systems used in en¬ 
ergy storage. NDCs are favored in various applications, 
such as hybrid devices and personal electronics due to 
their short load cycle and longer life. These properties 
have even greater advantage where flexibility and smaller 
sizes are desired. 

Since singe layer, hexagonal boron nitride (BN) can 
be grown on top of graphene,!i^Hi^ the use of BN as 
the solid dielectric nanoscale precision spacer between 
graphene layers is promising candidate for the fabrication 
of graphene based NDCs. On the other hand, graphene 
layers having perfect electron-hole symmetry can func¬ 
tion as metallic plates to store and release charge. Ear¬ 
lier, based on first-principles calculations, we predicted 
that an NDC consisting of a few insulating BN layers 
placed between two lattice-matched graphene layers can 
achieve high capacitance values.^i^ Furthermore, in that 
NDC model the capacitive behavior of the system di¬ 
verges from the classical Helmholtz model as the sepa¬ 
ration between the metallic plates gets smaller, due to 
the increase in the dielectric constant of BN thin films. 
This anomalous size-dependent increase of the capac¬ 
itance was also experimentally verified soon after our 
prediction.E^ Also it has been experimentally demon- 


(a) (b) 



FIG. 1: (a) PNDC formed by a zigzag BN stripe placed be¬ 
tween metallic graphene stripes, which display ID transla¬ 
tional periodicity along ^-direction. The unit cell is delineated 
by dashed lines, p and q are number of atoms in graphene 
and BN stripes in the unit cell, (b) Electronic band struc¬ 
ture of PNDC under zero bias voltage (or Ey=Q), where spin 
up and down states are shown with red(dark) and blue(light) 
lines, respectively, (c) Total charge density p(r) isosurfaces of 
PNDC. (d) {xz)- plane averaged electronic potential, V(y). 


strated that the permittivity of BN thin films increase 
significantly and leads to an increase in the capacitance 
value much higher than what is attained classically. 

In this study, we show that planar nanoscale dielec¬ 
tric capacitors (PNDCs) can be realized as a one-atom- 
thick, single layer honeycomb structure consisting of 
a BN stripe as dielectric spacer between two metallic 
graphene stripes. These laterally stacked stripes are lat- 
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FIG. 2: (a) Evolution of the energy bands of PNDC[4/4/4] under applied electric field Ey. Shifts of relevant spin bands are 
indicated by arrows. For Ey > 0.6 V /A, PNDC becomes a half-metal (HM). (b) A perspective view of the array of PNDC[4/4/4] 
periodically repeating along y-axis. (c) The isosurface of difference charge density, A(p), showing the charge separation, where 
the right graphene stripe is depleted from electrons, which are in turn deposited to left graphene stripe due to the shifts of 
bands under Ey. (d) (a; 2 )-plane averaged electronic potential V{y,Ey) exhibiting a potential difference of AV{y)=3V under 
Ey = l V/ A. 


tice matched. Furthermore, using these PNDCs, one can 
achieve high energy storage, as well as high gravimetric 
capacitance values, which are comparable to those of su¬ 
percapacitors. As an alternative to NDCs consisting of 
vertically stacked dielectric layers capped by two metal¬ 
lic graphene, PNDCs achieve charge separation between 
two separated, parallel metallic stripes embedded in the 
same atomic plane. 

It was shown experimentalljff^®^ that monolayer 
metallic graphene and dielectric BN can be grown com- 
mensurately with desired periodicity and with sharp 
boundaries. Thus one can construct not only single pla¬ 
nar capacitors but also various series and parallel com¬ 
binations of the former to achieve higher potential dif¬ 
ference or charge separation. Thus PNDCs predicted 
in the present study will make an important forward 
step towards the fabrication of atomically thi n circ uitry 
based on graphene/BN lateral heterostructure^^nm] to al¬ 
low microwave, as well as heat and pressure sensing ap¬ 
plications. 


II. METHOD 

Due to nanoscale sizes of the stripes, the stored energy 
through charge separation has to be treated from the 
first-principles. However, available first-principles meth¬ 
ods allow us to treat only one kind of uniform excess 
charge (negative or positive) in the same system at a 
time. The crucial aspect of this study is to maintain the 
charge separation of opposite polarity in two graphene 


stripes at both sides of BN by applying an external elec¬ 
tric field, Ey, in the plane of PNDC but perpendicular to 
the axis of the stripe. However, the electronic potential 
under the applied electric field makes a dip in the vac¬ 
uum spacing between PNDC stripes treated within the 
periodic boundary conditions using supercell geometry. 
As the strength of the electric field increases, this dip 
is further lowered as if a quantum well and allows the 
plane wave basis set to have states confined to the well 
as the solution of the Hamiltonian. Once the energies of 
these confined states are lowered below the Fermi level 
they start to be occupied by electrons. As a result, the 
electrons in the graphene layers are going to spill into the 
vacuum region under the external applied electric field. 
The spilling of the charge from PNDC to the vacuum is 
clearly erroneous and unrealistic. These artifacts of plane 
wave basis set become even more critical for wide vacuum 
spacing. On the other hand narrow vacuum spacing is 
not desired since it gives rise to significant coupling be¬ 
tween adjacent PNDC treated using the supercell geom¬ 
etry. These artifacts can be eliminated by using basis set 
consisting of orbitals centered only at the atomic sites, 
which does not allow the stat es co nfined to the poten¬ 
tial dip in the vacuum spacing.l^^l^ Thus, we carried out 
first-principles, spin-polarized calculations within density 
functional theory, where the eigenstates of Kohn-Sham 
Hamiltonian are expressed as linear combination of nu¬ 
merical atomic orbitals. The exchange-correlation po¬ 
tential was approximated by Perdew, Burke and Ernz- 
erhof functional.!^ A 200 Ryd mesh cut-off was chosen. 
Core states were replaced by norm-conserving, non-local 
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Trouiller-Martins pseudopotentialsP^ Atomic positions 
and lattice constants were optimized using the conju¬ 
gate gradient method by minimizing the total energy and 
atomic forces for each configuration. Dipole corrections^ 
were applied in order to remove spurious dipole interac¬ 
tions. All numerical calculations were performed using 
the SIESTA code.l^ 


III. RESULTS AND DISCUSSIONS 


Our model of PNDC described in Fig. 0^) is 
composed of lateral stacking of two narrow, zigzag 
graphene nanoribbons (stripes), ZGNR(p) and a zigzag 
BN nanoribbon, ZBNR(q) in between. Here p and q 
indicate the number of C and B-N atoms in their unit 
cells, respectively. Bare ZGNR(p) and ZBNR(q) are 
metallic. When saturated by hydrogen aton^ ZBNR(q) 
nanoribbon transforms into an insulator.i23 Similarly, 
when joined to graphene nanoribbons commensurately 
at both sides, the same ZBNR(q) can function as a di¬ 
electric spacer, while the graphene nanoribbons continue 
to be metal. These stripes can be set to desired widths 
by changing p and q values. While these three lattice 
matched graphene/BN/graphene composite stripes pre¬ 
serve their translational periodicity along x-axis and ac¬ 
quire ID character, the translational symmetry is bro¬ 
ken along y-direction. This model is treated by using 
periodic boundary conditions, where ID PNDGs are re¬ 
peated periodically along y- and z-axis keeping 10 A vac¬ 
uum spacing between them. The electronic structure 
of PNDG[p/q/p] structure under zero bias AU=0 (or 
Ey=0) is shown in Fig. [^b). It has spin polarized ground 
state with magnetic moment /i= 2.3 y,B- Graphene-like 
and BN-like bands are separated in the direct space and 
Fermi level is pinned at graphene sides while BN re¬ 
mains as an insulator. As seen in Fig. 0c), owing to 
the different ionicities of G, N, and B; a minute charge 
is transferred across the zigzag interface even for Ey=0. 
However, this does not lead to a capacitive behavior. In 
fact, the (xz)-plane averaged electronic potential V (y) in 
Fig. 0d) indicates a vanishing potential difference, i.e. 
AV (^=0 between graphene stripes. 

Normally, when connected to a real electric circuit, 
separation of charges of opposite polarity is attained be¬ 
tween graphene stripes of PNDC under the applied ex¬ 
ternal bias voltage AU, whereby the chemical potentials 
at metallic graphene stripes shift by eAU. In this non¬ 
equilibrium state an electric field is generated across the 
BN stripe. Here for the reason explained above, we sim¬ 
ulate this normal operation of PNDC by applying an 
electric field Ey along the y-direction with magnitude 
in the range of 0 V/A to 1 V/A. Notably, even if Ey=l 
v/A appears to be high, in earlier studies such high 
field values were attained.^ For each value of applied Ey, 
we first carried out calculations to optimize the atomic 
structure and lattice constants of PNDC by minimizing 


the total energy and atomic forces. Also, for each value 
of applied Ey we calculated the electronic energy bands 
of optimized structures as well as the net charge on the 
left graphene stripe, — Q in e/cell. This is obtained from 
the integral 


- Q = J Ap{r,EJy)dr (1) 

computed in the part of the PNDC unit cell correspond¬ 
ing to the left graphene stripe, where the boundary be¬ 
tween the left (right) graphene and BN stripe bisects C-B 
(N-C) bonds. Here the difference charge density corre¬ 
sponding to a given applied field Ey is calculated as 


Ap{r,Ey) = p{r,Ey) - p{r,Ey = 0) (2) 

where the total charge density of the PNDC under Ey=0 
is subtracted from the total charge density under the ap¬ 
plied Ey, and 


occ 

p{r, Ey) = -e ^ «'^.fc.<T(i’, F4)T„,fc,^(r, Ey). (3) 

n,k,(T 

Here, '^n,k,a- indicates the spin polarized electronic 
state. The same analysis was also performed for the right 
graphene stripe which is charged positively. It should be 
noted that the total charge density, p{r,Ey) integrated 
over the entire unit cell has to be zero, which has been 
verified also for each value of applied Ey. 

According to the atomic model presented in Fig. [^and 
Fig.0 the states of two graphene stripes have small cou¬ 
pling and hence their bands are split in momentum space. 
The evolution of the energy bands with applied electric 
field Ey is presented for 0.2 V/A< Ey < l.OV/A in 
Fig. 0a). Two unoccupied spin down bands (indicated 
with small arrows) shift in opposite directions, where the 
one associated with the left graphene stripe moves down 
and starts to become occupied. Concomitantly, a par¬ 
tially occupied spin up band associated with the right 
graphene stripe rises and hence becomes gradually de¬ 
populated. This situation explains how the charge and 
spin separations are achieved. Interestingly, for Ey ~ 0.6 
V/A the PNDC becomes a half-metal with p,=2.0 ps- 
Similar half-metallic behavior in charged graphene/BN 
composites were also reported before.^lEII The shorting 
of graphene stripes and hence the shorting of PNDC is 
hindered by the dielectric BN spacer between graphenes. 

In Fig. 0c) we present the difference charge density, 
Ap{y,Ey) of the PNDC under the field of Ey=\ V/A. 
The electron transfer from the right graphene stripe to 
the left stripe and the charge separation is clearly seen. 
In Fig. 0d), we show (a;z)-plane averaged electronic po¬ 
tential V{y). The potential at the left graphene stripe 
is lowered relative to the right graphene stripe, which 
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FIG. 3: (a) Variation of excess charge Q (e/cell); (b) stored 

energy Eg (eV/cell) and (c) the corresponding gravimetric 
capacitance C (F/g) for PNDC[4/4/4]. (d) Capacitance C 
values in (F/g) of the planar nanoscale dielectric capacitor 
PNDC[p/q/p] are calculated for p=4-12 and q=4-12 for spe¬ 
cific value of Ey for which C saturates. For each value of 
p, lower line connected by dots corresponds to capacitance 
values calculated through the expression, C = jlraEg^ 
while the upper line connected by crosses is computed from 
C' — QjmtEV. The calculated variation of C with tensile 
strain, e is shown by inset. 

is in compliance with the normal operation of PNDC. 
The electronic potential difference between two graphene 
stripes is calculated as AV = 3 V under Ey=l V/A. 

In Fig. [^a-b) we present the variation of calculated 
charge Q and stored energy Eg as a function of Ey. The 
energy stored on PNDC as a function of Ey, is obtained 
from the expression 

Eg{Ey)=ET{Ey)-ET{Ey=0), (4) 

where ET{Ey) is the total energy calculated for a given 
applied field Ey. Using the values of Q and Eg calculated 
for a given Ey, the capacitance per unit mass C can be 
obtained from 

C = Qy2mEg, (5) 

where m is the mass of PNDC per unit cell. The variation 
of the capacitance as a function of Ey is shown in Fig. 
j^c), which is high at small Ey (or bias voltage V) and 
gradually saturates to a constant value. 

To explore the effect of size of PNDC[p/q/p] we re¬ 
peated our analysis for different widths of graphene and 
BN stripes, namely for different p and q. Moreover, be¬ 
cause of uncertainties involved in the effective widths of 
metallic and dielectric stripes (hence in the boundaries of 


the integral given in equation]^, which affect the calcu¬ 
lated values of Q and C, here we present two alternative 
approaches in calculating C. In the first approach, C is 
calculated from the expression C — /2mEg, where 
Q and Eg correspond to Ey leading to the saturated 
values of C as shown in Fig. [^c). In the second ap¬ 
proach, we followed the definition of the capacitance as 
C = QlmAV. Results obtained from these two alter¬ 
native calculations are presented in Fig. [^d), where the 
capacitances of PNDCs are plotted for different q values. 
For both approaches, C decreases as q increases. Also 
C decreases as the width of graphene stripe p, increases. 
However, the difference between the values of capacitance 
calculated by using two alternative approaches, namely 
C — C increases for small p and q, but decreases and 
eventually diminish as p and q get larger. 

The capacitance values of PNDCs are found to de¬ 
pend strongly on the external conditions. For example, 
as shown by inset in Fig. id), C [4/4/4] changes by ^ 
30% under the tensile strain of e=10%, which can be 
detected easily and may offer technological applications. 
Similarly, when deformed under normal force or pressure, 
a PNDC can be deformed, whereby geometric parame¬ 
ters such as the effective width of BN, w is modified. 
Also, the energy barrier, <f> for electrons tunneling from 
left graphene to right graphene between two graphene 
stripes can be modified with local forces. Consequently, 
the capacitance value of PNDC is modified by local force 
or pressure. Capacitance value can also vary with the 
ambient temperature, since charge separation and hence 
Q, can be modified by the Fermi-Dirac distribution at 
finite temperature. Consequently, PNDCs can function 
as pressure and temperature sensors, when embedded in 
a soft materials like peptides or artificial tissues. It is 
noted that graphene/BN/graphene junctions of PNDC 
with small w can work also as tunneling diodes or res¬ 
onant tunneling device, where the tunneling current It 
can be monitored by mechanical deformation, bias volt¬ 
age or by temperature. Accordingly, these junctions also 
attain critical functionalities. Notably, all these features 
discussed above can be realized in diverse geometries, in 
particular in tubular forms. 

Next we consider a series combination of PNDC[p/q/p] 
which is composed of commensurate sequential stripes, 
graphene(p)/BN(q) /graphene(2p) /BN(q)/graphene(p) 
stacked laterally. Accordingly, the middle graphene 
stripe has the twice width of the graphene stripes at 
either side of series PNDC. Under the applied electric 
field the excess negative and positive charges were 
accumulated at the far left and right edges of the 
middle graphene(2p). The capacitance of this series 
combination was calculated as half of the capacitance 
of PNDC[p/q/p], namely Cgq = C[p/q/p] / 2. This 
demonstrates that to attain the desired Q or Eg or AU 
various combinations PNDC can be achieved on a single 
honeycomb structure by arranging metallic interconnects 
between parallel graphene stripes (which are graphene 
nanoribbons by themselves). 
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Finally, the PNDC structure revealed in this paper sug¬ 
gests us the construction of composite materials com¬ 
posed of the lateral stacking of a sequence of lattice 
matched zigzag ZGNR(p) / ZBNR(q) / ZGNR(p')... or 
armchair AGNR(p) / ABNR(q) / AGNR(p').. nanorib¬ 
bons embedded on a large single layer, 2D honeycomb 
structure. For p=4 and q=4 this structure can be 
viewed as line compound consisting of the adjacent nar¬ 
row stripes of graphene and BN, or (5-doping for p >> 4 
and q=4 or vise versa. In the composite materials with 
sequence p >> 4 and q » A the graphene stripes be¬ 
have as if multichannel ID metals on a 2D plane, which 
are separated by wide BN stripes. Depending on the val¬ 
ues of p and q, these composites can display a variety of 
electronic structures and thus herald the engineering of 
a large class of materials. 


IV. CONCLUSIONS 


In conclusion, we propose a planar nanoscale dielectric 
capacitor model which allow diverse parallel, series and 
mixed combinations to achieve desired values of charge 
separation, energy storage and potential difference val¬ 
ues, and can be fabricated in the same plane of micro or 
nanocircuits for crucial electronic applications. They are 
flexible and can also form in tubular geometry to serve as 
pressure, heat and strain sensors. Present results let us 
also infer interesting 2D composite materials constructed 
from commensurate and lateral stacking of graphene and 
BN nanoribbons with zigzag or armchair edges. 
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